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Fractal structures of normal and anomalous diffusion 
in nonlinear nonhyperbolic dynamical systems 

N. Korabefl and R. KlagesS 
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A paradigmatic nonhyperbolic dynamical system exhibiting deterministic diffusion is the smooth 
nonlinear climbing sine map. We find that this map generates fractal hierarchies of normal and 
anomalous diffusive regions as functions of the control parameter. The measure of these self-similar 
sets is positive, parameter-dependent, and in case of normal diffusion it shows a fractal diffusion coef- 
ficient. By using a Green-Kubo formula we link these fractal structures to the nonlinear microscopic 
dynamics in terms of fractal Takagi-like functions. 

PACS numbers: 05.45.-a, 05.60.-k, 05.10.-a, 05.45.Ac, 05.45.Df 



It is well-known that diffusion processes in nonlinear 
dynamical systems may be generated by microscopic de- 
terministic chaos in the equations of motion. This fact 
points to a somewhat deeper foundation of noncquilib- 
rium statistical mechanics than modeling diffusive trans- 
port by stochastic random walks. In order to understand 
deterministic diffusion, statistical mechanics was suitably 
combined with dynamical systems theory Jl|-|j|. Much 
was learned by analyzing simple models of deterministic 
transport such as one- and two-dimensional maps [^-0, 
chaotic billiards in external fields S, periodic Lorentz 
gases |^,|l0| , and certain differential equations [[ll| . How- 
ever, characteristics of deterministic diffusion were also 
observed experimentally, that is, in dissipative systems 
driven by periodic forces such as Josephson junctions in 
the presence of microwave radiation [i2| ], in superionic 
conductors [ fl3| , and in systems exhibiting charge-density 
waves [|l4j . The equations of motion of these systems are 
generally of the form of some nonlinear pendulum equa- 
tion. In the limiting case of strong dissipation, these 
differential equations were reduced to nonhyperbolic one- 
dimensional maps sharing certain symmetries |fl5|| . The 
so-called climbing sine map is a well-known example of 
this class of maps Due to its nonhyperbolicity, the 

map posseses a rich dynamics consisting of chaotic diffu- 
sive motion, ballistic dynamics, and localized orbits. Un- 
der parameter variation these different types of dynamics 
are highly intertwined resulting in complicated scenarios 
related to the appearance of periodic windows ||. On 
the other hand, for simple one-dimensional hyperbolic 
maps it was shown that the diffusion coefficient is typi- 
cally a fractal function of control parameters 16 r7j|. An 
analogous behavior was also detected for other transport 
coefficients pi], and in more complicated models [|s|T10||. 
However, up to now the fractality of transport coefficients 
could be assessed for hyperbolic systems only, whereas, 
to our knowledge, the fractal nature of classical transport 
coefficients in the broad class of nonhyperbolic systems 
was not discussed. 

Here we show that the nonhyperbolicity of the climbing 



sine map does not destroy these fractal characteristics of 
deterministic diffusive transport. On the contrary, fractal 
structures appear for normal diffusive parameters as well 
as for anomalous diffusive regions. We argue that higher- 
order memory effects are crucial to understand the origin 
of these fractal hierarchies in nonhyperbolic systems. By 
using a Green-Kubo formula for diffusion, the dynamical 
correlations are recovered in terms of fractal Takagi-like 
functions. We furthermore show that the distribution of 
periodic windows forms devil's staircase-like structures as 
a function of the parameter and that the complementary 
sets of chaotic dynamics have a positive measure in pa- 
rameter space that increases by increasing the parameter 
value. 

The climbing sine map we study is defined as 



M a (x n ) , M a (x) := x + asin(27ra) , (1) 



where a G R is a control parameter, and x n is the 

position of a point particle at discrete time n. Obviously, 
M a {x) possesses translation and reflection symmetry, 

M a (x + m) = M a {x) + m , M a (-x) = —M a (x) . (2) 

The periodicity of the map naturally splits the phase 
space into different cells (m, m + 1], m £ Z. We will 
be interested in parameters a > 0.732644 for which the 
extrema of the map exceed the boundaries of each cell 
for the first time indicating the onset of diffusive motion. 

The bifurcation diagram of the associated reduced map 
M a (x) :— M a (x) mod 1 , x := x mod 1, consists of in- 
finitely many periodic windows, see Fig. ^. Whenever 
there is a window the dynamics of Eq. (^) is either ballis- 
tic or localized || . Fig. [l] demonstrates that this scenario 
has a strong impact on the diffusion coefficient defined 
by D(a) := lim„^ 00 (a;^}/ (2n), where the brackets denote 
an ensemble average over moving particles. For localized 
dynamics orbits are confined within some finite interval 
in phase space implying subdiffusive behavior for which 
the diffusion coefficient vanishes, whereas for ballistic mo- 
tion particles propagate superdiffusively with the diffu- 
sion coefficient being proportional to n. Only for normal 
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diffusion D(a) is nonzero and finite. At the boundaries 
of each periodic window there is transient intermittent- 
like behavior eventually resulting in normal diffusion with 
D(a) r~j a^a) |^|||. Here we are interested in the com- 
plete parameter-dependent diffusion coefficient. For this 
purpose we compute D(a) from numerical simulations by 
using the Green-Kubo formula for maps P, p3||ll||l8| ] , 
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FIG. 1. Upper panel: bifurcation diagram for the climb- 
ing sine map. Lower panel: diffusion coefficient computed 
from simulations as a function of the control parameter a in 
comparison with the correlated random walk approximation 
D\ (a) (dots). The dashed vertical lines connect regions of 
anomalous diffusion, D(a) — > oo or D(a) —+ 0, with ballistic 
and localized dynamics in respective windows of the bifurca- 
tion diagram. 

D n (a) = (j a (i )J:(x)) - \{jl(x )) , (3) 

where the angular brackets denote an average over the in- 
variant density of the reduced map, (. . .} := J dxp(x) . . .. 
The jump velocity j a is defined by j a (x n ) := [x n +i] - 
[x n ] = [M a (x n )], where the square brackets denote the 
largest integer less than the argument. The sum J™ (x) := 
^^ =0 j a (ifc) gives the integer value of the displacement 
of a particle after n time steps that started at some ini- 
tial position x = xq, and we call it jump velocity func- 
tion. Eq. ([|) defines a time-dependent diffusion coef- 
ficient which, in case of normal diffusion, converges to 
D{a) = linin^oo D n (a). In our simulations we trun- 
cated J™ (x) after having obtained enough convergence 
for D(a), that is, after 20 time steps. The invariant den- 
sity was obtained by solving the continuity equation for 
p{x) with the histogram method of Ref. [EJ. 

The highly non-trivial behavior of the diffusion coeffi- 
cient in Fig. [j] can qualitatively be understood as follows: 
The Green-Kubo formula Eq. (j^) splits the dynamics into 
an inter-cell dynamics, in terms of integer jumps, and 
into an intra-cell dynamics, as represented by the invari- 
ant density. We first approximate the invariant density 
in Eq. (ffl to p(x) ~ 1 irrespective of the fact that it is 
a complicated function of x and a || . This approximate 
diffusion coefficient we denote with a superscript in Eq. 



(||), D^ip). The term for n = is well-known as the 
stochastic random walk approximation for maps, which 
excludes any higher-order correlations ||,||,[l7) • The gen- 
eralization -D*(a) , n > was called correlated random 
walk approximation |icj . We now use this systematic ex- 
pansion to analyze the diffusion coefficient of the climbing 
sine map in terms of higher-order correlations. 
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FIG. 2. (a) Sequence of correlated random walks D^(a) for 
n — 1, ■ ■ ■ , 10. Note the quick convergence for normal diffu- 
sive parameters. The dashed lines define the same periodic 
windows as in Figjl]. The inserts (b) and (c) contain blowups 
of D{ (a) in the initial region of (a). They show self-similar 
behavior on smaller and smaller scales. 

In Fig|| (a) we depict results for D n (a) at n = 
1, . . . , 10. One clearly observes convergence of this ap- 
proximation in parameter regions with normal diffusion. 
Indeed, a comparison of D\ (a) with D(a), as shown in 
Fig.|l|, demonstrates that there is qualitative agreement 
on large scales. On the other hand, for parameters cor- 
responding to ballistic motion the sequence of D^a) di- 
verges, in agreement with D(a) — > oo, whereas for local- 
ized dynamics it alternates between two solutions. This 
oscillation is reminiscent of the dynamical origin of lo- 
calization in terms of certain period-two orbits. That 
these solutions are non-zero is due to the fact that the 
invariant density was approximated. In regions of nor- 
mal diffusion this approximation nicely reproduces the 
irregularities in the diffusion coefficient. Even more im- 
portantly, the magnifications in Fig.|| give clear evidence 
for a self-similar structure of the diffusion coefficient. Our 
results thus show that dealing with correlated jumps only 
yields a qualitative understanding of normal and anoma- 
lous diffusion in the climbing sine map. 

We now further analyze the dynamical origin of these 
different structures. According to its definition, the time- 
dependent jump velocity function J™(x) fulfills the recur- 
sion relation 

JS(x)= Ja (x) + JS~ 1 (M a (x)). (4) 

J™(x) is getting extremely complicated after some time 
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steps, thus we introduce the more well-behaved function 




FIG. 3. Functions T™(x) for the climbing sine map rep- 
resenting the number of jumps particles starting at x have 
traveled when this number is integrated from zero to x, see 
Eqs. (Q)-(Q). (a) depicts diffusive dynamics at a — 1.2397 
(upper curve) and at a — 1.7427 (lower curve), (b) ballistic 
dynamics at a = 1.0, and (c) localized dynamics at a = 1.5. 
In (a) the limiting case for n — > oo is shown, in (b) and (c) it 
was n — 5, 6, 7. Note the divergence in (b) and the oscillation 
in (c). 

Integration of Eq. (ji|) thus yields the recursive functional 
equation 

t?(x) = t a (x) + -J— Tr\Mam - m ^ 

M a (x) 
with the integral term 

r M a (.x) 

m ■■= / dzg"{z)T^-\z) , (7) 

Jo 

where t a (x) := J dz j a {z), and g (z) is the second deriva- 
tive of the inverse function of M a (x) [19J. For piecewise 
linear hyperbolic maps I(x) simply disappears and the 
derivative in front of the second term reduces to the local 
slope of the map thus recovering ordinary de Rham-type 
equations ||[l7],[l8| . It is not known to us how to directly 
solve this generalized de Rham-equation for the climbing 
sine map, however, solutions can alternatively be con- 
structed from Eq. (ra) on the basis of simulations. Results 
are shown in Fig. pT For normal diffusive parameters the 
limit T a (x) = lim, woo T™(i) exists, and the respective 
curve is fractal over the whole unit interval somewhat re- 
sembling (generalized) fractal Takagi functions J3|,|l7],|l8| . 
However, in case of periodic windows T™(x) either di- 
verges due to ballistic flights, or it oscillates indicating 
localization. Interestingly, in these functions the corre- 
sponding attracting sets appear in form of smooth, non- 
fractal regions on fine scales, whereas the other regions 
appear to be fractal. 



The diffusion coefficient can now be formulated in 
terms of these fractal functions by integrating Eq. (||). 
For a G (0.732644, 1.742726] we get 

D(a) = 2 [T a (x 2 )p(x 2 ) - T a (xi)p(xi)] - D p a (a), (8) 

where Xi, i = 1,2, is defined by [M a (xi)] := 1, and 
Dq(cl) :— dxp(x). Our previous approximation 
Dl l (a) is recovered from this equation as a special case. 

The intimate relation between periodic windows and 
the irregular behavior of the diffusion coefficient moti- 
vates us to investigate the structure of the periodic win- 
dows in the climbing sine map in more detail. The ap- 
pearance of windows was analyzed quite extensively for 
non-diffusive unimodal maps |2Cf1 , whereas for diffusive 
maps on the line, apart from the preliminary studies of 
Refs. ||, nothing appears to be known. The windows are 
generated by certain periodic orbits, consequently there 
are infinitely many of them, and they arc believed to be 
dense in the parameter set Windows with ballistic 
dynamics are born through tangent bifurcations, further 
undergo Feigenbaum-type scenarios and eventually ter- 
minate at crisis points. Windows with localized orbits 
only occur at even periods. They start with tangent bi- 
furcations and exhibit a symmetry breaking at slope-type 
bifurcation points. 
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FIG. 4. Devil's staircase-like structure formed by the dis- 
tribution of periodic windows as a function of the control pa- 
rameter. N is the integrated number of period six-windows. 
The inset shows a blowup of the initial region. 



In order to analyze the structure of the regions of 
anomalous diffusion, we sum up the number of period 
six-windows as a function of the parameter, that is, the 
total number is increased by one for any parameter value 
at which a new period six-window appears. This sum 
forms a devil's staircase-like structure in parameter space 
indicating an underlying Cantor set-like distribution for 
the corresponding anomalous diffusive region, see Fig. ^. 
The (Lebesque) measure of periodic windows is obviously 
positive, hence this set must be a fat fractal pl| . Its 
self-similar structure can quantitatively be assessed by 
computing the so-called fatness exponent p2j. We are 
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furthermore interested in the parameter dependence of 
this fractal structure, therefore we divide the parameter 
line into subsets labeled by the integer value of the map 
maximum on the unit interval, [M a (x max )\ = j, j G Z. 
For j — 1,2,3 we obtain a fatness exponent of 0.45 with 
errors of 0.03, 0.04, and 0.05 for the different j. We men- 
tion that this value was conjectured to be universal and 
was also obtained for non-diffusive unimodal maps [ p2[ . 

We now study the measure of the windows as a func- 
tion of the parameter. For this purpose we computed 
all windows up to period six for the first subset, up to 
period 5 for j — 2, 3, and we summed up their mea- 
sures in the respective subsets. We find that the total 
measure decays exponentially as a function of j while os- 
cillating with odd and even values of j on a finer scale 
p3f . This oscillation can be traced back to windows gen- 
erated by localized dynamics, which only appear at even 
periods thus contributing only periodically to the total 
measure. However, different measures of 'ballistic' and 
'localized' windows decay with the same rate. We have 
furthermore computed the complementary measure Cj of 
diffusive dynamics in the jth subset of parameters. We 
find that Ci = 0.783, C 2 = 0.808, and C 3 = 0.932 with 
an error of ±0.002, so the measure of the diffusive re- 
gions is always non-zero and seems to approach one with 
increasing parameter values. 

We finally remark that the climbing sine map is of 
the same functional form as the respective nonlinear 
equation in the two-dimensional standard map, which 
is considered to be a standard model for many physical, 
Hamiltonian dynamical systems. Indeed, both models 
are motivated by the driven pendulum, both are strongly 
nonhyperbolic, and although the standard map is area- 
preserving it too exhibits a highly irregular parameter- 
dependent diffusion coefficient. Understanding the ori- 
gin of these irregularities was the subject of intensive re- 
search however, so far the complexity of the system 
prevented to reveal its possibly fractal nature. A suit- 
ably adapted version of our approach to nonhyperbolic 
diffusive dynamics may enable to make some progress in 
this direction. 
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